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Computer: DEC Station 5000/260; other machines with FORTRAN 77 and suf- 
ficient capacity should be able to run this program 

Installation: Institut for Theoretical Physics, University of Frankfurt, Frank- 
furt/Main 

Operating system: tested under UNIX but does not depend on the particular 
operating system 

Programing language used: FORTRAN 77 

High speed storage required: 5 Mbytes 

Number of bits per byte: 8 

Number of lines in combined test deck: 19296 

CPC subprograms used: Jetset 7.4 

Keywords: polarized nucleon-nucleon scattering, high energy physics, Monte 
Carlo simulation 

Nature of physical problem: This program can be used to simulate polarised 
nucleon - nucleon collisions at high energies. Spins of colliding particles are 
taken into account. The program allows to calculate cross sections for various 
processes. 

Method of solution: The existing Monte Carlo program Pythia 5.7 has been 
modified to incorporate spin effects. The program incorporates all features 
of Pythia. 

Restrictions on the complexity of the problem: Spins of colliding nucleons must 
be parallel to the collision axis. 

Typical running time: ^0.1 sec CPU per event on DEC Station 5000/260 

References: Reader has to refer to Pythia manual for additional details. 
Further informations can be obtained from: Prof. A. Schafer, Institut fiir 
Theoretische Physik, Universitat Frankfurt a. M., Robert Mayer-Strasse 10, 
D-60054 Frankfurt am Main, Germany 
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LONG WRITE-UP 



1 Introduction 

The detailed investigation of hadronic spin effects, most notably in deep- 
inelastic lepton-nucleon scattering, has proven to allow for extremly sensi- 
tive tests of QCD. One of the most interesting issues is the possible con- 
tribution of the axial-vector anomaly. This anomaly plays a central role in 
modern quantum-field-theory and spin effects might offer the only possibil- 
ity to actually detect it. To estimate its contribution one has to know the 
distribution function for polarised gluons A(yf(x, Q^)- It is argued that the 
total spin carried by gluons could be large, because as{Q'^)^g{Q'^), with 
^giQ"^) = lo ^9(^1 Q'^)i is a renormalisation group invariant, such that if 
AglQ"^) is a few percent at the constituent quark level it could be very size- 
able at large Q^. In deep- inelastic scattering it is impossible to distinguish 
between the regular quark contribution and the gluon contribution via the 
anomaly. There is actually no conceptually sound way to separate them. In 
exclusive or semi-inclusive lepton-hadron scattering channels it is in princi- 
ple possible to discern a gluonic and a quark contribution, but the by far 
most direct and conceptually clean way are polarised proton-proton colli- 
sions. RHIC [H offers the option to study polarised proton collisions up to 
y/s = 300 GeV which would be ideal for this purpose. To analyse this option 
in detail, to plan the experiments, and to determine the reachable precission 
a detailed Monte Carlo program is needed, which we provide with this pub- 
lication. This program can handel, however, only longitudinal polarisation. 
All polarisation effects in proton-proton reactions share the same problem. 
As only a small fraction of the partons is polarised the signal to background 
ratios are typically of the order percent, such that a detailed simulation is 
necessary to extract quantitative results. 

Because a state-of-art Monte Carlo program is usually the result of many 
years of development we modified an existing, widly used program called 
P YHTIA |l|] . We have been interested only in those parts of the program that 
correspond to nucleon-nucleon collisions. Consequently only these routines 
were changed. The resulting polarised version of Pythia is named Sphinx, 
an acronym for Simulator of Polarised Hadronic YNteractp^jions. 

It is interesting to note that Sphinx is able to function in unpolarised 
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mode as well. In that case it is functionally equivalent to Pythia. This 
means that results will be the same but because some data structures inside 
the program have been modified the program requires more space. 

The main problem was to introduce the changes in a way compatible 
with the basic structure of Pythia, e. g. the way in which Pythia evaluates 
the cross sections or simulates the effects of higher order QCD corrections. 
Pythia is based on the parton model. The main assumption is that the 
nucleon may be represented by an incoherent mixture of quarks and gluons 
carrying portions of nucleon momentum. Their momentum distributions 
are described by parton distributions. In the parton model the nucleon is 
composed of quarks and antiquarks of various fiavours and of gluons. In the 
polarised case each parton comes in one of two helicities: righthanded and 
lefthanded. Therefore the number of parton components has to be doubled, 
corresponding to the two degrees of freedom for longitudinal polarisation. 
The most natural idea was to generalize the notion of fiavour to account 
for the polarisation. So instead of d-quark we have lefthanded d-quark and 
righthanded d-quark etc. Because in Pythia gluons are treated as a special 
"fiavour" the scheme applies to them as well. As far as the general structure 
of the program is concerned the only complication of the polarised code is 
that we have more "fiavours" and this is rather easy to implement. 

Polarisation is followed in Sphinx only up to the hard partonic inter- 
action, i. e. the hadronic cross section and the higher order corrections in 
the initial state shower are evaluated spindependently whereas all final state 
interactions as fragmentation and decays are treated spinindependently. No 
theoretical model for polarised fragmentation exists. It is, however, generally 
accepted that it should not depend on longitudinal polarisation. For trans- 
verse polarisation the situation would, however, be different. One knows 
experimentally that substantial effects exist, e. g. a strong coupling between 
transverse spin and transverse momentum. 

The polarisation effects have been implemented in the following parts of 
the program: 

1. Parton distributions. At the moment we supply six sets of polarised 
parton distributions: Altarelli-Stirling two parametrisations from 
a Ross- Roberts article and three sets from Gehrmann&Stirling [^]. 
First order Altarelli-Parisi evolution is taken into account in all 
cases. 
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2. Hard processes. The processes currently implemented in the polarised 
mode are summarized in Table |I|. Other processes may be used as well 
but the results will be averaged over polarisations. The cross sections 
for all polarised processes are given in the literature [0 , and we checked 
them. 

3. Initial state showering. The polarised case was not discussed in this 
context before. We obtained all necessary formulae. 

4. Documentation. The event listing has been modified such that it dis- 
plays polarisation information about the interacting particles up to the 
hard partonic interaction. 

Hereafter we do not try to explain the structure of Pythia 5.7 and 
we indicate only the most significant modifications which lead to Sphinx. 
Readers unfamiliar with Pythia should first consult the appropriate manual 

i- 

2 From Pythia to Sphinx: modifications 

Our program is an extension of Pythia version 5.7. The Pythia 5.7 pro- 
gram heavily relies on Jetset 7.4 subroutines ||^. These subroutines are 
mainly used to simulate processes that take place after the hard partonic 
interaction, i.e. fragmentation and decays, or provide spinindependent ma- 
nipulations as for example Lorentz-transformations. Because polarisation 
effects in the final state are neglected in Sphinx, the Jetset program could 
be left untouched. Only the set-up of an event listing which is also a task of 
Jetset in a Pythia run has been taken over from Sphinx in the polarised 
mode. The event listing in Sphinx provides then in addition information 
about the polarisation flow. For that purpose two new subroutines have 
been added which are modifications of the corresponding Jetset routines 
(see below). By this means Sphinx provides a suitable interface to Jet- 
set such that the programs may be linked without any modifications in the 
latter. There is only one exception: the Jetset subroutine LUGIVE, which 
returns the values of all Jetset and Pythia common block variables cannot 
be longer used, because in Sphinx the dimensions of some Pythia arrays 
have been enlarged so that the formats no longer match. 
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Table 1: List of processes implemented in the polarised mode. 



I SUB 


Process 


Comment 


1 




quark-antiquark annihilation into virtual 


2 




annihilation into charged vector boson 


11 


QiQi QiQi 


(anti-) quark - (anti-) quark scattering; anni- 
hilation diagram not included 


12 


QiQi qkQk 


annihilation process 


13 


QiQi 99 


annihilation into gluon pair 


14 


QiQi 91 


annihilation into gluon and prompt 7 


15 


QiQi gZ^ 


annihilation into gluon and Z^ 


16 


QiQi gW^ 


annihilation into gluon and 


18 


qiqi 11 


annihilation into 7 pair 


19 


q,q, -jZ^ 


annihilation into 7 and Z° 


20 


(Mi ";^r= 


aiiiiihilatioii into and 11'^ 


28 


qi9 (ii9 


(anti-)quark - gluon scattering 


29 


qi9 Qii 


prompt 7 production in (anti-)quark - gluon 

scattering 


30 


qi9 QiZ° 


Z^ production in (anti-)quark - gluon 
scattering 


31 


qi9 QjW^ 


production in (anti-)quark - gluon 

scattering 


53 


99 QkQk 


gluon fusion 


68 


99 99 


gluon - gluon scattering 
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In many places inside Pythia partonic data are stored in arrays indexed 
by flavour. Inclusion of helicities means doubling the number of "flavours" 
as described in Section |I[ For that reason we had to enlarge the arrays. The 
unpolarised mode has been preserved throughout the program and works 
exactly as the original code. All modifications are implemented in a treelike 
structure with several branching points where the user may choose between 
the polarised and unpolarised modes. The default is always unpolarised. An 
interesting point is that one can combine polarised and unpolarised treat- 
ments for different processes and for various stages of the same process. This 
has to be done with care and one has to be aware of this not always being 
justified from the physical point of view. The branching points were imple- 
mented by means of the IF-ELSE structure which directs the program flow 
according to a specific switch. Each functional part (usually a subroutine) 
has its own local variable called IPOL. This variable controls the mode and 
its value depends on some parameter^ MSTP(x) and NSUB(x) to be defined 
by the user in the main program. For example, it is possible to run the code 
with polarised partonic cross section and polarised parton distributions but 
unpolarised initial state showering etc. This solution provides flexibility but 
has to be used carefully. The major disadvantage of IF-ELSE constructions 
is that parts of the code are multiplied. This makes servicing of the code 
more cumbersome because one has to introduce the same change in different 
places at once. In addition the resulting code became quite lengthy. 

Sphinx is as Pythia a "slave-system", i. e. it consists only of callable 
subroutines where two of them (PYINIT and PYEVNT) have to be called by 
the user to perform the event generation and a few other, e. g. DPLIST, 
PYSTAT, etc., could be called to obtain further event information. The rest 
of the subroutines is of internal use. Therefore the user has to supply a main 
program where all relevant parameters and switches have to be specified and 
these subroutines have to be called. Because this structure is the same as 
in Pythia the reader is referred to again for details about the general 
Pythia parameters. In this article only the new parameters in Sphinx 
are discussed. The modifications made in the individual subroutines are 
described in Section |[ In Table |l| we list all new parameters introduced to 
control the polarised mode. Finally in Section ^examples for a main program 
and the result of the corresponding test runs are presented. 

iSee Table 
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Table 2: Parameters controlling the polarised mode 



Parameter 


Description 


Default 


MSTP(171) 


=0 
=1 
=2 


beam polarisation 
unpolarised 

polarisation in +z direction 
polarisation in —z direction 





MSTP(172) 


=0 
=1 
=2 


target polarisation 
unpolarised 

polarisation in +z direction 
polarisation in —z direction 





MSTP(175) 


use of polarised parton distribution in po- 
larised initial state shower (MSTP(176)=l) 

=0: unpolarised distribution; for testing only, do 
not use! 

=1: polarised distribution 


1 


MSTP(176) 


initial state showering mode 
=0: unpolarised 
=1: polarised 





MSTP(177) 


set of polarised parton distributions Aq and 
Ag used; in addition one has to specify un- 
polarised set as in standard Pythia 

=0: Ag = and A^f = (no polarisation) 

=1: fake polarisation, built up from unpolarised 
distribution according to Ag = "^^xoo^^'' ^ 

=2: Altarelli-Stirling parametrization [Q; data 
file altsti.dat required 

=3: Ross- Roberts parametrization |^ set d 

=4: Ross-Roberts parametrization set a; data 
file rosroa.dat required 

=5: Gehrmann- Stirling parametrization set a; 
data file partons.dat required 

=6: Gehrmann- Stirling parametrization set b; 
data file partons .dat required 

=7: Gehrmann- Stirling parametrization [^] set c; 
data file partons . dat required 

=8: fake polarisation, for testing only, do not use! 
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Table 2: Parameters controlling the polarised mode, continued 



Parameter 



Description 



MSTP(178) 



percentage of fake polarisation for 
MSTP(177)=1 



MSTP(180) 



mode selection (master switch) 
=0: unpolarised mode; this value overrides all 

other polarisation switches 
=1: polarised mode 



NSUB(ISUB) 



mode for subprocess ISUB 
=0: unpolarised treatment 
=1: polarised treatment 



NSEL 



menue of polarised processes 



=1 


ISUB = 


11 , 12 , 13 , 28 , 53 , 68 switched on 


=10 


ISUB = 


14,18, 29 switched on 


=11 


ISUB = 


1 switched on 


=12 


ISUB = 


2 switched on 


=13 


ISUB = 


15,30 switched on 


=14 


ISUB = 


16,31 switched on 
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Table 3: Internal variables storing polarisation information 



Variable 


Description 


Com. Block 


MINTOll) 


=2 


beam helicity 
unpolarised 
positive helicity 
negative helicity 


PYINTl 


MINT(312) 


=0 
=1 
=2 


target helicity 
unpolarised 
positive heUcity 
negative helicity 


PYINTl 


MINT(313) 


=0 
=1 
=2 


helicity of shower initiator on beam side 
unpolarised 
positive helicity 
negative helicity 


PYINTl 


MINT(314) 


=0 
=1 
=2 


helicity of shower initiator on target side 
unpolarised 
positive helicity 
negative helicity 


PYINTl 


MINT(315) 


=0 

=1 
=2 


helicity of hard interacting parton on beam 
side 

unpolarised 
positive helicity 
negative helicity 


PYINTl 


MINT(316) 


=0 

=1 
=2 


helicity of hard interacting parton on target 
side 

unpolarised 
positive helicity 
negative helicity 


PYINTl 


MSTP(179) 


switch off polarisation temporarely in PYSIGH 

and PYSTFU resp. 
=0: no action 
=1: switch off polarisation 


PYPARS 
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Table 3: Internal variables storing polarisation information, continued 



Variable 



Description 



Com. Block 



ISIGH(1000,6) 

ISIGH(I,1) 
ISIGH(I,2) 
ISIGH(I,3) 
ISIGH(I,4) 
ISIGH(I,5) 
ISIGH(I,6) 



hard scattering information of Ith line 
particle code of Ith hne on beam side 
particle code of Ith hne on target side 
colour flow 

helicity of Ith line on beam side 
helicity of Ith line on target side 
not used 



PYINT3 



KD(I) 



polarisation/hehcity of Ith line 
no polarisation/helicity 
positive polarisation/helicity 
negative polarisation/helicity 



DPYPOL 



XSFX(2, -40:40, 0:2) 



XSFX(JT,KFL,0) 
XSFX(JT,KFL,1) 

XSFX(JT,KFL,2) 



X times parton distribution for given 
X and of flavour KFL = -40:40 
and helicity KFLD = 0:2 on beam side 
(JT=1) and target side (JT=2) resp. 
unpolarised 
positive helicity 
iiegati\'e hc4ic-it}' 



PYINT3 
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3 Common Blocks and Subroutines 



In the following the modified and the new subroutines that have been specif- 
ically created for Sphinx are described in more detail. The general tasks of 
the subroutines themselves as well as the unchanged parameters and variables 
are not explained, because they are the same as in Pythia. The reader is 
asked again if occasion arises to consult |I| to obtain the needed information. 
We restrict our explanation to the new aspects in Sphinx. The purpose of 
the modifications is indicated and the new parameters, switches, and internal 
variables are listed. Meaning and possible values of the new parameters are 
given in Table |^. To incorporate polarisation the following common blocks 
have been enlarged and replace the corresponding Pythia common blocks 
or are added: 

• C0MM0N/PYINT3/XSFX(2,-40:40,0:2) ,ISIG(1000,6) , SIGH (1000) 

• C0MM0N/PYSUBS/MSEL,NSEL,MSUB(200) ,NSUB(200) , KFIN (2 , -40 : 40) , 
CKIN(200) 

• C0MM0N/DPYP0L/KD(4000) 

Information about the new internal variables and enlarged arrays can be 
found in Table 0. In addition it is shown how the local polarisation switch 
IPOL is built up in the different subroutines. Only the polarised case (IPOL=l) 
will be discussed, because in the unpolarised case (IPOL=0) each subroutine 
works exactly as the corresponding Pythia subroutine. The not mentioned 
subroutines of Sphinx are the same as in Pythia. 

MAIN PROGRAM 

Purpose: to set up the polarised event generation. The variables which 
have to been set are listed in Table |l|. 

Remarks: Examples of a main program are given in Section ^. 
SUBROUTINE PYINIT (FRAME , BEAM , TARGET , WIN) 

Purpose: to display Sphinx header; to check partially the availability of 
the desired polarisation scenario, i. e. to control that the master switch for 
polarisation MSTP(180) is set properly, the selected partonic subprocesses 
can be treated polarised and to control and compose the polarisation menue 
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via NSEL; to call DPLIST instead of LULIST (see below). 
New parameters: MSTP(180), NSEL, NSUB(ISUB) 
Internal Polarisation switch: IP0L=MSTP(180) 

Remcirks: If a not allowed scenario has been chosen, the programs stops 
with an appropriate error message. 

SUBROUTINE PYEVNT 

Purpose: to start polarised event generation; to call DPEDIT instead of 

LUEDIT (see below). 

New parameters: MSTP(180) 

Internal polarisation switch: IP0L=MSTP(180) 

SUBROUTINE PYINKI (CHFRAM , CHBEAM , CHTARG , WIN) 

Purpose: to check availability of the desired hadronic polarisation scenario, 
i. e. to control that the selected hadron can be treated polarised and to 
verify that the polarisation is longitudinal; to store the polarisation of beam 
and target for the event listing in KD(1) and KD(2) and for internal use in 
MINT(311) and MINT (312). 

New parameters: KD(I), MSTP(171), MSTP(172), MSTP(180) 
New internal variables: MINT (311), MINT (312) 
Internal polarisation switch: IP0L=MSTP(180) 

Remarks: At present only nucleons and hyperons and their antiparticles 
can be treated polarised. If a not allowed scenario has been chosen, the 
programs stops with an appropriate error message. 

SUBROUTINE PYRAND 

Purpose: to adapt PYRAND to the new environment - all relevant arrays 
which have been enlarged or added to the common blocks in other subrou- 
tines are modified here as well; to extend event shape selection to incorporate 
helicities; to store helicities of the partons entering the hard interaction ac- 
cording to 

• MINT(313): helicity of the beam parton for use in the initial state 
showering subroutine; 

• MINT (314): helicity of the target parton for use in the initial state 
showering subroutine; 
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• MINT(315): helicity of the beam parton; 

• MINT (316): helicity of the target parton. 

New parameters: MSTP(180), NSUB(ISUB) 

New internal variables: MINT(313), MINT(314), MINT(315), MINT(316) 
Internal polarisation switch: IP0L=MSTP(180) xNSUB(ISUB) 
Remarks: Note that MINT (313) =MINT (315) and MINT (314) =MINT (316) 
but the values of MINT (313) and MINT (314) are changed later by the ini- 
tial state shower in PYSSPA. 

SUBROUTINE PYSCAT 

Purpose: to adopt PYSCAT to the new environment (see PYRAND); to store 

helicities of the partons entering the hard interaction; to fill lines 1, 2 and 5, 

6 in the event listing with polarisation information (see below). 

New parameters: KD(I), MSTP(180), NSUB(ISUB) 

New internal variables: MINT(315), MINT(316) 

Internal polarisation switch: IP0L=MSTP(180) xNSUB(ISUB) 

SUBROUTINE PYSSPA (IPUl , IPU2) 

Purpose: to perform polarised initial state showering, hehcity dependent 
GLAP evolution equations are used in the backward evolution algorithmus; 
to enlarge all relevant array in an appropriate manner to incorporate polari- 
sation; to check proper selection of the polarised initial state shower scenario, 
i. e. to control that the MSTP(175) and MSTP(176) are set correctly; to store 
the hehcities of initial state shower initiators (MINT (313), MINT (314)). 
New parameters: KD(I), MSTP(171), MSTP(172), MSTP(175, MSTP(176), 
MSTP(180), NSUB(ISUB) 

New internal variables: MINT (313), MINT (314), MINT (315), MINT (316) 
Internal polarisation switch: IPOL=MSTP ( 180) x NSUB ( ISUB) x MSTP (176) 
Remarks: At the present stage only QCD shower can be treated polarised, 
QED showering has to be done in the unpolarised manner. The combination 
MSTP ( 175) =0 and MSTP ( 176) =1 allows to simulate polarised showering with 
the use of unpolarised parton distributions. This option is just for testing 
and should not be selected by the user! If MSTP (175) or MSTP (176) are set 
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improperly, the programs stops with an appropriate error message. The in- 
ternal variables MINT (313) and MINT (314) are changed to their final values 
in this subroutine. 

SUBROUTINE PYMULT(MMUL) 

Purpose: to switch off polarisation in PYSIGH (set MSTP(179)=1 temporally) 
when called from PYMULT even in a polarised run, because multiple interaction 
cannot be treated polarised at the moment. 
New parameters: MSTP(179) 

SUBROUTINE PYREMN(IPU1 , IPU2) 

Purpose: to adopt PYREMN to the new environment (see PYRAND); to fill lines 
3, 4 in the event listing with polarisation information (see below). 
New parameters: KD(I) 

SUBROUTINE PYSIGH 

Purpose: to evaluate the helicity dependent hadronic cross sections by con- 
volution of the helicity dependent parton distributions with the helicity de- 
pendent partonic cross sections; to supply the subroutine with the helicity 
dependent partonic cross sections. 

New parameters: MSTP(171), MSTP(172), MSTP(179), MSTP(180), 
NSUB(ISUB) 

Internal Polarisation switch: 
IP0L=MSTP(180) xNSUB(ISUB) X (1-MSTP(179) ) 

Remarks: PYSIGH will always run in the unpolarised mode when it is called 
by PYMULT which sets MSTP(179)=1 in PYSIGH temporally. Evaluating the 

spindependent hadronic cross sections one has to notice that the hadrons are 
specified according to the spin, whereas the partons are labelled by their helic- 
ities. The spin is defined relative to the collision axis and the beam is assumed 
to move in the positive direction. For that reason the helicities at the target 
side are opposite to the polarisations. Hence the target labels are reversed in 
the convolution in comparison to the beam labels. The parton distributions 
are passed from PYSTFU to PYSIGH through the array XPQ(KFL,KFLD) (see 
below) and stored in the array XSFX(N,KFL,KFLD), where KFLD denotes the 
helicity. ISIG(N, I) contains the information about the Nth-line in the event 
listing. The new entries ISIG(N,4) and ISIG(N,5) store the helicities of the 
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partons at the beam and target side respectively. ISIG(N,6) is reserved but 
not used at the moment. 

SUBROUTINE PYSTFU(KF,X,Q2,XPQ) 

Purpose: to evaluate the helicity dependent parton distributions for given 
flavour (KF), x (X), and (Q2) according to the selected parametrisations. 
New parameters: MSTP(177), MSTP(178), MSTP(179), MSTP(180), 
NSUB(ISUB) 

Internal Polarisation switch: IP0L=MSTP(180) x (1-MSTP(179) ) 
Remarks: Call of PYSTFU returns x times the parton distribution functions 
for given flavour, x, and for both helicities and an averaged (unpolarised) 
value. The values are stored in the array XPQ(KFL,KFLD) which has been 
enlarged from XPQ(-25:25) to XPQ(-25:25,0:2). Row XPQ(KFL,0) con- 
tains the unpolarised distributions, XPQ(KFL,1) distributions corresponding 
to the positive helicity (relative to the hadron) and XPQ(KFL,2) distributions 
for the negative helicity. The parton distributions are selected by switches 
described earlier (see Table p. The polarised distributions g-t, g± are con- 
structed from the unpolarised ones g, g (selected by old Pythia switches) 
and polarised parts Ag and lS.g selected by MSTP(177) (see Table |I]). Four 
subroutines have been added to calculate the polarised distributions. These 
are ALTSTlQ, ROSROA, ROSROD, and GEHSTI0. The subroutines require data files 
altsti.dat, rosroa.dat, and partons.dat. These files must be visible to 
FORTRAN open statement and therefore they have to be placed in appropriate 
directory. The files are supplied with the program. The parton distribution 
for fixed helicity are reconstructed from the polarised and unpolarised distri- 
butions according to g-t = | (g ± Ag). When the polarised and unpolarised 
parts are combined together the program performs the unitarity check - if 
the resulting total distribution becomes negative for one helicity it is put to 
zero and the corresponding result for the other helicity is set to the value 
of the unpolarised part. Only polarised parametrisations for protons are 
implemented. Neutron parametrisations are obtained from them by isospin 
symmetry, the parametrisations for hyperons are constructed by naive SU(3). 

■^This subroutine has been written by G. Altarelli and J. Stirling. Used with permision 
from the authors. 

^ The three parametrisation of Gehrmann&Stirhng contained in GEHSTI are brandnew 
and not thouroughly tested yet. 



16 



Charge conjugations is used to describe the corresponding antiparticles. 

SUBROUTINE ALTSTI (X , Q2 , UPV , DNV , SEA , STR , CHM , BOT , TOP , GLU) 

Purpose: to return x times the polarised parton distributions evahiated at 
given X (X) and (Q2) according to the parametrisation of Altarelh&Stirhng 
[0]. UPV denotes the valence distribution of up-quarks xAu^^^{x, Q"^), DNV for 
down-quarks xA(i^'^^(a;, Q^)- SEA signifies the sea distribution xAq^'^'^{x,Q'^). 
STR, CHM, BOT, and TOP label the distributions for the strange-, charm-, 
bottom-, and top-quark xAg(a;,Q^), q = s,c,b,t respectively. Finally GLU 
marks the gluon distribution xAglxjQ"^). 

Remarks: ALTSTI requires the data file altsti .dat which has to be placed 
in an appropriate directory, altsti.dat is supplied with this program. 

SUBROUTINE R0SR0D(X,Q2,XPDF) 

Purpose: to return x times the polarised parton distributions evaluated at 
given X (X) and (Q2) according to the parametrisation of Ross&Roberts 
set d 0. XPDF(-6 : 6) contains xAq{x, Q^) for q = t,b, c, s, u, d, g, d, u, s, c, b, t 
in this order. 

SUBROUTINE R0SR0A(X,Q2,XPDF) 

Purpose: to return x times the polarised parton distributions evaluated at 
given X (X) and (Q2) according to the parametrisation of Ross&Roberts 
set a [0. The contents of XPDF(-6:6) is explained above. 
Remarks: ROSROA requires the data file rosroa.dat which has to be placed 
in an appropriate directory, rosroa.dat is supplied with this program. 

SUBROUTINE GEHSTI (IGFLAG , X , Q2 , XDDPR) 

Purpose: to return x times the polarised parton distributions evaluated at 
given X (X) and (Q2) according to the parametrisation of Gehrmann&Stir- 
lingset a (IGFLAG=0), set b (IGFLAG=1) or set c (IGFLAG=2). [§. XDDPR(-6:6) 
contains xAg(x, Q"^) for q = i,b, c, s, u, d, g, d, u, s, c, b, t in this order. 
Remarks: GEHSTI requires the program package parton. f, consisting of 
the subroutines polpar, parini, getpar, q21ow, q2high, xlow, and xhigh 
and the data file partons.dat written by Gehrmann and Stirling and used 
here with permission of the authors. For further informations about this 
package see 0. partons.dat has to be placed in an appropriate directory 
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and is supplied with this program. The package part on. f is contained in 
this program. 

SUBROUTINE DPLIST(MLIST) 

Purpose: to display the polarisations of the particle in the event hsting. 
New parameters: KD(I) 

Remarks: DPLIST is a modification of the Jetset subroutine LULIST. It is 
changed to display the polarisation in the final listing. The sign displayed 
just behind the particle code denotes polarisation with respect to the z-axis. 
When the sign is missing the particle has been treated as unpolarised. The 
information is taken from the vector KD(I) and transformed accordingly to 
('0','l'j'2')~^(' '/+';'"')• The following format is chosen (polarisation for the 
colliding hadrons, helicity for the partons resp.): 

• = ' ': no polarisation/helicity 

• = '+': positive polarisation/hehcity 

• = '— ': negative polarization/helicity 
SUBROUTINE DPEDIT(MEDIT) 

Purpose: to compress the vector KD(I), containing the polarisation infor- 
mation, properly. 
New parameters: KD(I) 

Remarks: DPEDIT is a modification of the Jetset subroutine LUEDIT. 

4 Examples 

In the following we give two examples of a main program for a simulation 
with Sphinx and show the corresponding results. We considered longitudi- 
nal polarised proton-proton scattering in the CMS at = 200 GeV and 
selected the process qg — > qg for the partonic interaction. In the first ex- 
ample both beam and target are polarised in +2;-direction, in the second 
example the target spin is reversed. With regard to the event listings is has 
to be mentioned that the displayed format differs from the real because we 
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removed a few columns such that it fits in this text. In addition the event 
hsting is cut after the hard interaction (denoted by • • •), i. e. behind fine 8. 
In the omitted part there is no polarisation information and it has the same 
format as the original Pythia listing. 

The information about the polarisation flow is displayed as explained 
above right behind the flavour code KF. The flrst row contains the information 
about the beam particle. In addition to the Pythia labels the sign '+' 
behind the flavour code for the proton KF=2212 denotes the polarisation in 
positive 2;-direction. Accordingly the sign '+' (' — ') in the second row signifles 
the polarisation of the target in positive (negative) ^-direction in the flrst 
(second) example. The third and fourth line represents the initiators of the 
initial state shower. In both examples they are gluons with positive helicities 
for both the beam and the target side. During the initial state shower the 
beam side parton becomes an s with positive helicity, whereas the target 
side parton remains a positive helicity gluon. These partons undergo the 
hard interaction the resulting partons of which are displayed in the lines 
seven and eight. In the following the flnal state interaction takes place, 
i. e. the outgoing partons fragment, the unstable produced hadrons decay, 
etc. as long as only stable particles exist. In the final state polarisation 
is not traced and consequently there is no polarisation informations about 
these lines provided. This part of the listing is then again the same as the 
corresponding Pythia listing. 

4.1 The Main Programs 

4.1.1 First Example — parallel polarisation 

C Example of a Main Program for event generating 

C in longitudinal polarised proton-proton-scattering 



C, 
C, 
C, 
C, 
C, 
C, 
C, 
C, 



E.g. : polarised p(+)p(+) -scattering in CMS 
at sqrt(s)=200 GeV 



This program has to be linked with 

the programs SPHINX and JETSET7.4, 

the data files ALTSTI.DAT and ROSROA.DAT, 

and the CERN-Libraries . 
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C COMMON BLOCKS of SPHINX for event generation 

C0MM0N/LUDAT1/MSTU(200) ,PARU(200) ,MSTJ(200) ,PARJ(200) 
COMMON/LUJETS/N , K (4000 , 5) , P (4000 , 5) , V (4000 , 5) 
C0MM0N/PYSUBS/MSEL,NSEL,MSUB(200) ,NSUB(200) , 
& KFIN(2, -40:40) ,CKIN(200) 

C0MM0N/PYPARS/MSTP(200) ,PARP(200) ,MSTI(200) , PARI (200) 
C0MM0N/PYINT5/NGEN (0 : 200 , 3) , XSEC (0 : 200 , 3) 

C Polarisation Set-Up 

C=============================================================: 

C Polarised simulation 

MSTP(180)=1 

C Beam positive polarised 

MSTP(171)=1 

C Target positive polarised 

MSTP(172)=1 

C Polarised parton distributions a la Altarelli&Stirling 

MSTP(177)=2 

C Polarised Initial State Shower 

MSTP(176)=1 
C with polarised parton distributions 

MSTP(175)=1 

C=============================================================: 

C=============================================================: 

C Event Set-Up 

C=============================================================: 
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C Choice of parton processes ''a la carte'' 

MSEL=0 

C Choice of polarised parton processes ''a la carte'' 

NSEL=0 

C Choice of process: qg — > qg 

MSUB(28)=1 

C Process 28 polarised 

NSUB(28)=1 

C Number of generated events 

NEVENT=1000 

C kinematical cuts 

C___ P_T-cut (minimum) 

CKIN(3)=3. 
C P_T-cut (maximum) 

CKIN(4)=25. 

C========================================================== 

c========================================================== 

C Start of event generation 

C========================================================== 

C Initialisation 

CALL PYINITCCMS' , 'p' , 'p' ,200.) 

C Loop over events 

DO 100 I=1,NEVENT 

C Event generation 

CALL PYEVNT 
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List the first event with spininf ormation 
IF(I.EQ.l) CALL DPLIST(l) 



100 CONTINUE 

C Print cross section and histogram 

CALL PYSTAT(l) 

END 

4.1.2 Second Example — antiparallel polarisation 

The second example is constructed by replacing MSTP(172)=1 by MSTP(172)=2 
in the first example, i. e. by switching the spin of the target. 

4.2 The Event Listings 
4.2.1 First example — peirallel polarisation 

SPHINX 

** Last date of change: 6 Aug 1994 ** 

The Lund Monte Carlo - PYTHIA version 5.7 
** Last date of change: 3 Apr 1992 ** 

The Lund Monte Carlo - JETSET version 7.4 
** Last date of change: 10 Mar 1992 ** 

1********** PYINIT: initialization of PYTHIA routines ********** 



I I 

I PYTHIA will be initialized for a p on p collider I 

I at 200.000 GeV center-of-mass energy I 

I I 
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*PYMAXI : summary of differential cross-section mcLximum search * 



I ISUB Subprocess name 
I 



I I 
I Maximum value I 
I I 



28 f + g -> f + g 

96 Semihard QCD 2 -> 2 



I 
I 
I 
I 



3.9475E+00 
1 . 7380E+02 



************** PYINIT: initialization completed *************** 



Event listing (summary) 



I 


particle/ jet KF orig 


P- 


.X 


P- 


-y 


P- 


.z 


E 


m 


1 


!p+! 


2212+ 








000 





000 


99 


996 


100.000 





938 


2 


!p+! 


2212+ 








000 





000 


-99 


996 


100.000 





938 


3 


!g! 


21+ 


1 





130 


-0 


241 


33 


098 


33.099 





000 


4 




21+ 


2 





151 


-0 


503 


-3 


866 


3.902 





000 


5 


!s~! 


-3+ 


3 


-0 


838 


-0 


101 


27 


030 


27 . 043 





000 


6 


!g! 


21+ 


4 


-0 


036 


2 


481 


-0 


350 


2.505 





000 


7 


!s~! 


-3 





2 


361 


2 


098 


26 


049 


26 . 240 





199 


8 


!g! 


21 





-3 


235 





281 





631 


3.308 





000 



1** PYSTAT: Statistics on Number of Events and Cross-sections *** 
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I Subprocess I Number of points I Sigma I 
I III 
I 1 1 (mb) I 

I III 
I N:o Type I Generated Tried I I 

I III 



I III 
I All included subprocesses I 1000 5855 I 6.695E-01 I 

I28f + g->f + g I 1000 5855 I 6.695E-01 I 

I III 



** Fraction of events that fail fragmentation cuts = 0.00000 *** 



4.2.2 Second example - antiparallel polarisation 

SPHINX 

** Last date of change: 6 Aug 1994 ** 

The Lund Monte Carlo - PYTHIA version 5.7 
** Last date of change: 3 Apr 1992 ** 

The Lund Monte Carlo - JETSET version 7.4 
** Last date of change: 10 Mar 1992 ** 

1********** PYINIT: initialization of PYTHIA routines ********** 



I I 

I PYTHIA will be initialized for a p on p collider I 

I at 200.000 GeV center-of-mass energy I 

I I 



*PYMAXI : summary of differential cross-section mciximum search * 
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I ISUB Subprocess name 
I 



I I 
I Maximum value I 
I I 



28 f + g -> f + g 

96 Semihard QCD 2 -> 2 



I 
I 
I 
I 



3 . 9772E+00 
1 . 7380E+02 



************** PYINIT: initialization completed *************** 



Event listing (summary) 



I 


particle/ jet KF orig 


P- 


.X 


P- 


-y 


P- 


.z 


E 


m 


1 


!p+! 


2212+ 








000 





000 


99 


996 


100.000 





938 


2 


!p+! 


2212- 








000 





000 


-99 


996 


100.000 





938 


3 


!g! 


21+ 


1 





784 





462 


33 


089 


33.102 





000 


4 




21+ 


2 


-0 


502 





119 


-33 


887 


33.891 





000 


5 


!s~! 


-3+ 


3 


-0 


240 





737 


21 


851 


21.865 





000 


6 


!g! 


21+ 


4 


-5 


178 


-4 


192 


4 


698 


8.152 





000 


7 


!s~! 


-3 





-4 


508 





563 


25 


310 


25.715 





199 


8 


!g! 


21 





-0 


910 


-4 


018 


1 


239 


4.302 





000 



1** PYSTAT: Statistics on Number of Events and Cross-sections *** 



Subprocess 



I I 
I Number of points I 
I I 



Sigma 



25 



I 
I 



I 
I 



I 
I 



(mb) 



I 
I 
I 
I 



I N:o Type 
I 



I Generated 
I 



Tried I 
I 



I I 



I I 

5870 I 6.868E-01 I 

5870 I 6.868E-01 I 

I I 




1000 
1000 



** Fraction of events that fail fragmentation cuts = 0.00000 *** 
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